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O I Abstract 



A simple inflationary model based on loop quantum cosmology is considered. 
Within this framework, we show that inflation does not necessarily erase the infor- 
mation prior to its onset, but that such information may leave its imprint in the 



;_i ' energy-spectrum of the gravitational-waves generated at these earliest of times. 

a , 

It is usually assumed that inflation, by its characteristics, among them a tremendous 
increase in the scale of the universe taking place in an extremely short time, will effectively 
remove any kind of information or properties already existing prior to the onset of this 
phenomenon (cf., e.g., |1] and references therein). In this work we use a simple example. 
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from the literature on loop quantum cosmology (LQC) [2]- [12], to show that this is not 
necessarily the casqj. In particular, from the model herein analysed, although the very first 
quantum stages in LQC are followed by inflation, still, traces of that pre-infiationary phase 
survive, leaving their imprint in today's energy-spectrum of the stochastic background of 
gravitational waves. 

Before proceeding into more technical matters, let us specify several elements that 
distinguish our framework. To calculate the power spectrum of the gravitational waves, 
we apply a method based on the Bogoliubov coefficients and the differential equations 
obeyed by them, a method first devised by L. Parker [H] and subsequently used in other 
papers [151 HSl El [IB]- Furthermore, our model includes several stages of expansion: a 
pre-infiationary era, followed by infiation and then the standard stage, bringing us to 
the present universe, whose dark energy component is, for simplicity, represented by a 
cosmological constant. In addition, we avoid using the sudden transition approximation 
[T8] and, instead, with the help of Parker's differential equations, consider the transitions 
between the different stages as continuous. 

In the model that we are considering, the potential U = —, appearing in the equa- 
tions for the gravitational waves (a designates the scale of the universe and the prime its 
derivative with respect to conformal time), shows, in the pre-inflationary phase, a specific 
feature, not present in the usual model based on the classical equation for the infiaton field. 
This feature is a small peak in the graph of U with respect to conformal time. This small 
barrier gives rise to a different set of initial conditions, for those gravitons whose frequen- 
cies satisfy k"^ < U . This, in turn, will be reflected in the shape of the energy-spectrum in 
the low-frequency range, for frequencies approximately below 10~^^ rad/s. For gravitons 
with frequencies above this, corresponding to energies above f/, the spectrum will not be 
changed, when compared to the classical case. 



We begin with a description of our background model, closely following [TT] and 
In LQC, a scalar fleld (f) with potential V{(f)), in a flat Friedmann- Robert son- Walker back- 
ground, is described by the Hamiltonian [11] 

n^ = aV(0) + ^d.^ipl, (1) 

where p^ = djj (p is the momentum canonically conjugate to (p and dj^i{a) = Di{q)/a? gives 
the eigenvalues of the geometrical density operator in loop quantization; j and / are the 
so-called ambiguity parameters, j being a half integer quantum number greater than one 
and where I determines the behavior of the density operator on small scales (compared to 
a fundamental scale a^,, introduced below) and can take any value between zero and one. 
Di{q) is a quantum correction factor for the density in the semiclassical regime [201 12] • 
Throughout the paper we take j as a free parameter. On the other hand, we flx / = 3/4. 



^Recently, a few papers addressing this or other related issues and frameworks, have appeared in the 
Hterature [131 



The expression for D-ii^i^q) is then given by 



D„M = (^) q^'\7[{q + ir/'-\q-l\''/']-llq[{q + iy"-sgn{q-l)\q-l\''/']}\ 



(2) 



where the variable q = [a/a^Y and a* = ■\/{'^j)/^lpi is a fundamental length scale arising 
in loop quantization; 7 is the Barbero-Immirzi parameter, for which calculations of the 
black-hole entropy give 7 ~ 0.274 [201 E]- 

The LQC corrections will appear in the background equations of our model through the 
function D{q) and its derivative t) (q). These effects are particularly relevant for q < 1, 
when the scale factor a < a^,. For g ^ 1 we enter the classical regime, characterized by 
D{q) — > 1 and D (q) — )■ 0. At this point the semiclassical regime reduces to the classical 
case. 

We assume a universe described by a flat FRW metric and going through three stages 
of expansion: 

(a) The pre-infiationary and the inflationary stages: To model it, we use the chaotic 
infiaton potential 

V{<P) = Imy. (3) 

The background equations are the Friedmann equation, H =a /a, 



rr2 _ 8^,2 



Vi^) + ^ 0^ 



and the modified Klein-Gordon equation for the infiaton, obtained from T-Ld 



DV^i^) = 0, 




(4) 



(5) 



where we take derivatives with respect to comoving time. The classical equations are 
obtained putting D = 1 and D= 0. 

(b) Transition between inflation and the radiation dominated universe: We de- 
fine the end of the inflationary era, beginning of the transition era, as the point 
where a becomes negative. During the transition, the energy of the infiaton will be 
transferred to the radiation fluid. We model this transfer through the action of a 
frictional term characterized by a decay constant F^. The equations above are now 
slightly changed. The Friedmann equation becomes [21] 



H' 
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Vi<P)+n^ +PR 



(6) 



and the equation for the scalar field is 



+3- (j) +V^{(P) = -T4 



(7) 



The equation for the radiation fluid is (we use F^ ~ 2 x 10 ^): 

Pn=-4-pn + r^(j) . (8) 

a 

At the beginning of the transition era we assume pr = 0. We end this period when 
when pi?a'^ = const. = PmOq. 

(c) Standard stage: This brings us to the present universe, whose dark energy compo- 
nent is taken, for simplicity, to be described by a cosmological constant. We, then, 
integrate the equations till the present time, with the content of the universe charac- 
terized by a radiation field p/j, dark matter+baryonic matter, with equation of state 
p = 0, and a cosmological constant representing the dark energy. The Friedmann 
equation is now 



H^ SttG 






Pro [ —J + Paid y— ) + Pdeo 



3 
the subscript denoting the values at the present time, a{to) = oq. 



(9) 



Above, we wrote our expressions using comoving time (as this is the most appropriate 
variable when we have the fast increase in the scale factor during inflation). However, the 
equations for the gravitational waves take a much simpler form when we use conformal 
time r, defined by writing the FRW metric in the form 

ds"^ = s2(r) {-dr^ + [5ij + h,j{r, x)] dx'dx^} . (10) 



The tensor perturbations hij can be expanded, in the usual manner, in terms of plane 

d^k 



waves: 



/d k 
,, ,.,,, , , rrr [«p(^, k)e,,(k,p)e^''-e(r, k) + h.c] 



11 



where x denotes the spatial coordinates, k = |k| = 2ti/\ = ua is the comoving wave 
number; the index p runs over the two polarizations of the gravitational waves and eij is 
the polarization tensor, a^ the annihilation operator and ^ the mode function. In conformal 
time, the mode function obeys the equation 

e + {k'-a"/a)^ = 0, (12) 

the derivatives being with respect to r. In what follows, we call gravitational wave potential 
the expression 

a" 
U = ^. (13) 

a 

This potential is shown in figure 1, for the semiclassical case (cf. eq. (I3])-(|5])). The LQC 
main feature appears as an additional small barrier at very early times, prior to the onset 
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Figure 1: Gravitational wave potential \U\, as it appears in the LQC model herein em- 
ployed. 



of the inflationary stage. This small barrier gives rise to different initial conditions at the 
end of the pre-inflationary era, for those gravitons whose frequencies /c^ < U. This, in 
turn, will be reflected in the shape of the energy-spectrum in the low-frequency range, as 
will be seen. 

More generically, for appropriate forms of a" {t) / a{T) we may have exponentially grow- 
ing solutions, with the gravitational field pumping energy into the gravitational waves. 
From the point of view of quantum mechanics, this corresponds to graviton produc- 
tion. In fact, equation ( 1T2|) closely mimics a Schrodinger equation with potential barrier 
U{t) = a" {t) / a{T) (with the time coordinate r instead of the spatial coordinate). It is an 
important point to realize that, whenever k"^ < U{t), the conditions exist for a significant 
graviton production, the equation for the gravitational waves being then of the form of a 
parametric oscillator, while, when k"^ ^ U{t), the equation for the gravitational waves is 
that of a harmonic oscillator, with no gravitons being produced. 

We now express the creation and annihilation operators in terms of the initial creation 
and annihilation operators v4jj(k) and Ap{]i.), through the Bogoliubov coefficients a{k,T) 
and/3(A;,r) [H], [22]: 



(r,k) = a(r,A;)Ap(k) + /3*(r,fc)At(k), 



(14) 



a and f3 satisfying the relation 

|a|'-|/3r = l. (15) 

These coefficients, which in the sudden transition approximation (cf. [HI [23]) are calculated 
by requiring the mode functions and their derivatives to be continuous across the transition. 



obey the following couple of differential equations 



a' = \\a + I3e^'^^-^^^A - (16) 

ZK '- -'a 



i r „ on^i^ ^ M a" 



and 

/3' = -^[/3 + «e"2^'=(— ')]^, (17) 

which, with a change to the variables X{k^ r) and Y{k^ r) 

a = ]-{X + Y)e'^^^-^^^ (18) 

/3 = ^(X-F)e-'^("-"^), (19) 

become 

X' = -ikY (20) 

completing our system of equations. These variables obey the constraint equation ( 1T^ . 
which takes the form 

XrYr + XiYi = l, (22) 

the subscripts r and z denoting the real and imaginary parts. In terms of the new variables, 

|/3P is 

m' = \ [{X^ - Y,f + (X, - Y^f] . (23) 

During the LQC dominated phase, we have, as mentioned, an extra barrier (cf. Fig. 
1). We verified numerically that the appropriate equation of state is p^ ~ p^, in which 
case an analytical solution of equations (1201) and fl^T]) can be found and is of the form 

x = ^[-Ue) + iY,{e)], (24) 

^ = \/^{-^o(^) + 2^^i(^) + ' ['^o(^) - 2^^i(^)]} , (25) 

V 0(7 

where J„ and Yn are the Bessel functions of the first and second kind, respectively. In 
this case, |/3p has a more complicated form, but it also satisfies the limits |/3p — >• 0, when 
k — )■ 00, and |/3p — >■ 00, when /c — )■ 0. Due to the smallness of the extra barrier, only 
the low frequency gravitons will be affected. For the higher frequency modes, only the 
inflationary barrier will be important. 

Initial conditions for the scale factor ai and for (pi need to be set, with (pi being obtained 
from (/)i with the help of the uncertainty principle (see [TT] ) 

10"^ / a,- ^ 



p/'^ \a^ 
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The scale factor is Oj > Ipi, otherwise the semiclassical equations would not be valid, 
due to the discretization of the geometry of space introduced by LQG. We therefore take 

We choose parameters and initial conditions such that we get enough inflation (number 
of e- folds ~ 70) and satisfy the large-scale CMB anisotropies requirements [1]: m^ ~ 
10~^?np and (pmax ^ Snip. We still have the ambiguity parameter j, as a free parameter. 

In the inflationary era we used, for convenience, Planck units and worked with the 
adimensional time /i = m^t. Finally, we used the same values for those parameters which 
are common to both models, m^ = 10~^, F^ = 2 x 10~^ and <pi/m^ = 2 for the semiclassical 
and 0i = 3.3 for the classical cases. 

It is now possible to integrate the full set of equations, including the equations for X 
and Y, fl2Up and f l2ip . The number of gravitons of gravitons created is, as is well known, 
given [T31|Tni[IB] by |/3/ina«P- Taking into account that the density of states is u'^du/{2n'^c^) 
and that each graviton contributes with two polarizations 2/la;, then, from the definition of 
the energy density dE = P[u)du, we have the following expression for the power-spectrum 

p{uy. 
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(27) 



We choose to present our results in terms of the dimensionless relative logarithmic energy- 
spectrum of the gravitational waves, at present time Tq : 
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Pcrit being the value of the present time critical density and pg^ the gravitational wave 
energy density 

Pg^= [ P{co)duj. (29) 




Figure 2: The energy-spectrum of the gravitational waves, for the LQC model employed 
in this paper, showing the sharp decrease in the low-frequency part of the spectrum. 

Our results for the energy-spectrum of the gravitational waves are plotted in figure 2. 
For purposes of illustration, we defined values for the parameters of the classical model 



such that the results would approximately comply with the total amount of anisotropy 
measured by COBE which, in terms of Q, is 

n < 1.37 X 10"^°, (30) 

for frequencies corresponding to the present horizon size Uhor, where we use H = 71 kms~^ 
Mpc~^, giving Uhor — 10~^^ rarfs"^. We immediately notice that the gradient of the low- 
frequency part of the spectrum is much steeper than in the classical case, a change that 
is due to the extra barrier in Figure 1. This shows that physical features present in the 
pre-inflationary eras are not erased out by inflation, as is commonly assumed, but may 
leave their imprint in the spectrum of the gravitational wavecl, which was the main point 
of our brief note. We leave to a future paper a more thorough investigation of this topic. 
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